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ABSTRACT 

The locations and amplitudes of three acoustic peaks and two dips in the Boomerang, 
MAXIMA and DASI measurements of the cosmic microwave background (CMB) 
anisotropy power spectra as well as their statistical confidence levels are deter- 
mined in a model-independent way. It is shown that the Boomerang-2001 data 
(Netterfield et al. 2001) fixes the location and amplitude of the first acoustic peak 
at more than 3a confidence. The next two peaks and dips are determined at a con- 
fidence level above la but below 2a. The locations and amplitudes of the first three 
peaks and two dips are ip^ = 212 ± 17, Ap^ = 5426 ± 1218 fiK^, la, = 413 ± 50, 
Ad^ = 1960 ± 503 fiK^, ip^ = 544 ± 56, Ap^ = 2266 ± 607 fiK^, Id^ = 746 ± 89, 
Ad^ = 1605 ± 650 ^iK^, ip.^ = 843 ± 35, Ap^ = 2077 ± 876 fiK"^ respectively {la errors 
include statistical and systematic errors). The MAXIMA and DASI experiments give 
similar values for the extrema which they determine. For MAXIMA these are the 1st 
and 3rd peaks, for DASI the 1st and 2nd peaks and the 1st dip. Moreover, the locations 
and amplitudes of the extrema determined from the combined data of all experiments 
are quite close to the corresponding values extracted from the Boomerang data alone. 

In order to use these data in a fast search for cosmological parameters an accurate 
analytic approximation to calculate CMB peak and dip positions and amplitudes in 
mixed dark matter models with cosmological constant and curvature is derived and 
tested. 

The determined cosmological parameters from the CMB acoustic extrema data show 
good agreement with other determinations, especially with the baryon content as de- 
duced from standard nucleosynthesis constraints (Buries et al. 2001). These data sup- 
plemented by constraints from direct measurements of some cosmological parameters 
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and data on large scale structure (LSS) lead to a best-fit model which agrees with 
practically all the used experimental data within la. The best-fit parameters are: 
= 0.64101^, Urn = O.Setg:?^, = 0.047lg:°^3, = 1.0lg:f?, h = 0.65l°:p and 
Tc = O.ISIq^j. The best-fit values of Qi, and T/S are close to zero, their la upper 
limits are 0.17 and 1.7 respectively. 

Subject headings: cosmology: microwave background anisotropics - acoustic peaks - 
cosmological parameters 

1. Introduction 

The new data on the cosmic microwave background (CMB) temperature anisotropy obtained 
in the Boomerang (de Bernardis et al. 2000; Netterfield et al. 2001), Maxima I (Hanany et al. 2000; 
Lee et al. 2001) and DASI (Halverson et al. 2001) experiments provide relatively accurate mea- 
surements of the CMB anisotropies up to £ ~ 1000. Boomerang is a long duration balloon (LDB) 
flight around the South Pole, MAXIMA is a balloon flight from Palestine, Texas, and DASI is an 
interferometer experiment. The mutual agreement of such divers experiments within statistical 
uncertainties is very reassuring. 

After a correction of the first results from Boomerang (de Bernardis et al. 2000) by Netter- 
field et al. (2001), these measurements are in astounding agreement with the simplest flat adiabatic 
purely scalar model of structure formation. The best fit cosmological parameters obtained coin- 
cide with other, completely independent determinations, like e.g. the baryon density parameter 
predicted by nucleosynthesis (Buries et al. 2001). 

CMB anisotropies can be calculated within linear perturbation theory in a multi-component 
universe. These calculations are very well established and allow accurate predictions of the CMB 
power spectrum for a given model of initial perturbations and given cosmological parameters. 
All the calculations are linear and very well controlled. Publicly available codes, e.g. CMBfast 
(Seljak & Zaldarriaga 1996) provide 1% accurate results for a given model within two minutes of 
CPU time on an ordinary PC. Due to these advantages, CMB temperature fluctuation data are 
extremely valuable for testing theoretical models of structure formation and for the determination 
of cosmological parameters. 

Nevertheless, the efficiency of parameter determination using codes like CMBfast to compute 
the temperature anisotropy spectrum for each model has several problems: 1) The complete set of 
observational data of the current state and the early history of the Universe is described by models 
with at least six parameters. The implementation of CMBfast-like codes into search procedures for 
best fits in high dimensional parameter spaces consumes too much CPU time even for the most ad- 
vanced computers due to the necessity to carry out numerical integration of the Einstein- Boltzmann 
system of equations which describe the evolution of temperature and density perturbations of each 
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component through the decouphng epoch. 2) The CMB power spectrum alone has several more or 
less exact degeneracies in parameter space (see e.g. Efstathiou & Bond (1999)) which can only be 
reduced substantially or removed completely if other data sets, e.g. galaxy clustering data, cor- 
responding to different scales and redshifts, arc combined with CMB measurements. The results 
and especially the error bars which are obtained from search procedures using different classes of 
cosmological observations with different quality and different statistical properties are difficult to 
interpret. 

Several groups have overcome the first problem by computing a grid of CMB anisotropy spectra 
in the space of models and interpolating between them to obtain the spectra for intermediate values 
of the parameters (see Tegmark et al. (2001); Lange et al. (2001); Balbi et al. (2000); de Bernardis 
et al. (2001); Wang et al. (2001) and references therein). 

Here we propose an alternative method: The CMB angular power spectrum obtained by COBE 
(Smoot et al.l992), Boomerang, MAXIMA-1 and DASI has well defined statistical and systematic 
errors in the range of scales from quadrupole up to the spherical harmonic £ ~ 1000 and the 
present data can be represented by a few dozen uncorrelated measurements. Practically the same 
information is contained in a few characteristics such as the amplitude and inclination of the power 
spectrum at COBE scale and the amplitudes and locations of the observed acoustic peaks and 
dips. Indeed, it was shown (see (Hu & Sugiyama 1995; Hu &; White 1996; Efstathiou & Bond 1999; 
Doran & Lilley 2001) and references therein) that the sensitivity of the locations and amplitudes 
of the extrcma, especially of the peaks, to cosmological parameters is the same as their sensitivity 
to the flat band powers, the C/'s, from presently available data. This is not surprising: the cxtrema 
have an obvious physical interpretation and have the largest weight among data points on CMB 
power spectrum as a result of a minimal ratio of error to value. Hu et al. (2001) have shown that 
most of the information from the Boomerang and MAXIMA-1 data sets can be compressed into 
four observables: amplitude and location of 1-st peak, and amplitudes of 2-nd and 3-rd peak. 

The first three acoustic peaks and the two dips indicated by the above mentioned experiments 
and the COBE large scale data can be presented by not more than 12 experimental points. If we 

use the approach by Bunn & White (1997) for the four year COBE data and the data on acoustic 
peaks, we have 7 experimental values to compare with theoretical models. Each of them can be 
calculated by analytical or semi-analytical methods. This enables us to study present CMB data 
in a very fast search procedure for multicomponent models. 

Even though, in principle, a measurement of the entire Ce spectrum contains of course more 
informations than the position and amplitudes of its peaks and dips, with present errorbars, this 
additional informations seems to be quite modest. 

The goal of this paper is to use these main characteristics of the CMB power spectrum to 
determine cosmological parameters. To do this we have to accomplish the following steps: 1) 
to locate the positions and amplitudes of three peaks and two dips as well as determining their 
error bars from experimental data, 2) to derive accurate analytical approximations to calculate 
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these positions and amplitudes and test them by fuU numerical calculations. We also derive an 
accurate and fast semi-analytical method to normalize the power spectrum to the 4-year COBE 
data. Such analytical approximations have been derived in the past for the matter power spec- 
trum (Eisenstein & Hu 1998; Eisenstein & Hu 1999; 1) and for the Sachs-Wolfe part of the CMB 
anisotropy spectrum (Kofman & Starobinsky 1985; Apuncvych &: Novosyadlyj 2000). Here we 
derive an analogous approximation for the acoustic part of the CMB anisotropy spectrum by im- 
proving an approximation proposed by Efstathiou &; Bond (1999). 

The outline of the paper is as follows. In Section 2 we determine the locations and amplitudes 
of the 1st, 2nd and 3rd acoustic peak as well as 1st and 2nd dip and their confidence levels using the 
published data on the CMB angular power spectrum from Boomerang (Netterfield et al. 2001), 
MAXIMA-1 (Lee et al. 2001) and DASI (Halverson et al., 2001). Analytical approximations for 
the positions and amplitudes of the acoustic peaks and dips are described in Section 3. A new 
method for an accurate and fast COBE normalization is also presented in this section. Details are 
given in two Appendices. Our search procedure to determine cosmological parameters along with 
the discussion of the results is presented in Section 4. In Section 5. we draw conclusions. 

2. Peaks and dips in the CMB power spectrum: experimental data 

We have to determine the locations and amplitudes of acoustic peaks and dips as well as their 
uncertainties in the data of the angular power spectrum of CMB temperature fluctuations obtained 
in the Boomerang (Netterfield et al. 2001), MAXIMA-1 (Lee et al. 2001) and DASI (Halverson 
et al. 2001). We carry out a model-independent analysis of the peaks and dips in the power spectra 
for each experiment separately, as well as using all data points jointly. 

2.1. Boomerang-2001 

A model-independent determination of peak and dip locations and amplitudes in the Boomerang 
data (Netterfield et al. 2001) has been carried out recently by dc Bcrnardis et al. (2001). Our ap- 
proach is based on a conceptually somewhat different method, especially in the determination of 
statistical errors. 

At first, mainly for comparison of our results with de Bernardis et al. (2001), we fit the peaks 
in the Boomerang CMB power spectrum by curves of second order (parabolas) as shown in Fig. 1. 
The six experimental points {N(,xp = 6) in the range 100 < £ < 350, which trace the first acoustic 
peak, are well approximated by a three parameter curve {Npar = 3). Hence, the number of degrees 
of freedom Nf = N^xp — Npar for the determination of these parameters is Nf = 3. The best-fit 
parabola has Xmm ~ extremum, located at £p-^ = 212 and Ap^ = 5426 nK'^, is the best-fit 

location and amplitude of the first acoustic peak. We estimate the statistical error in the following 
way. Varying the 3 parameters of the fitting curve so that ~ Xmin < 3.53 the maxima of the 
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parabolas define the la range of positions and amplitudes of the first acoustic peak in the plane 
{i, £{£ + l)Q/27r). The boundary of this region determines the statistical la errors for the location 
and amplitude of the first acoustic peak. We obtain 

In Fig. 1 the la, 2a (the boundary of the region with — Xmin — ^-O^) and 3a {x^ — Xmin — 14.2) 
contours are shown in the plane (£, l{i + l)(7^/27r). All contours for the first acoustic peak are closed. 
This shows that the BGomcrang-2001 data prove the existence of a first peak at a confidence level 
higher than 3a. The values Ax^ = 3.53, 8.02 and 14.2, given by the incomplete Gamma function, 
Q{Nf, Ax'^/2) = 1 - 0.683, 1 - 0.954 and 1 - 0.9973, correspond to 68.3%, 94.5% and 99.73% 
confidence levels respectively for a Gaussian likelihood oi Nf = 3 degrees of freedom. These levels 
which depend on Nf and thus on the number of independent data points (which we just took at 
face value from Netterfield et al. (2001)) define the regions within which the maxima of parabolas 
leading to the data points lie with a probability > p, where p = 0.68, 0.954 and 0.9973 for 1-, 2- 
and 3-0" contours respectively. The same method for the second peak using the eight experimental 
points in the range 400 < i < 750, hence A^^ = 5, gives 

ip2 = 541]']j^Q2, Ap^ = 2225^227 ■ 
For the third peak (750 < £ < 1000, 6 experimental points, Nf = 3) we obtain 

£p, = 843^11, Ap, = 2077^^^^ fiK\ 

For the second and third peaks, the 2- and 3-a contours axe open as shown in Fig. 1. This 
means that the Boomerang experiment indicates the 2nd and 3rd acoustic peaks at a confidence 
level higher than 1- but lower than 2-a. This is in disagreement with the result obtained in 
(de Bernardis ct al. 2001). Formally the disagreement consists in the fact that de Bernardis et 
al. (2001) set Nf = 2 for all peaks and dips (see paragraph 3.1.2 of their paper) leading to different 
values of Ax^ for the 1-, 2- and 3-a contours. They argue that there are two free parameters, namely 
the height and the position of the peak. The 'philosophy' of the two approaches is somewhat 
different: While our contours limit the probability that the given data is measured if the correct 
theoretical curve has the peak position and amplitude inside the contour, in their approach the 
contours limit the probability that the given best fit parabola leads to data with peak position and 
amplitude inside the contour. In other words, while they compare a given parabola to the best 
fit curve, we compare it to the data. In that sense we think that the closed 2-a contours of de 
Bernardis et al.(2001) do not prove the existence of the secondary peaks at the 2-a level. Of course 
our approach has a problem as well: It relies on the data points being independent. If they are 
not, the number of degrees of freedom should be reduced. 

The same procedure can be applied for the amplitudes and positions of the two dips between 
the peaks. 
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There is also a slight logical problem in the results presented in Fig. 2 of de Bcrnardis et 
al.(2001): If 2nd and 3rd peaks are established at 2cr C.L., then the 2nd dip should be determined 
at the same C.L.; but even the la contour for the position of the second dip is not closed. Also 
the position of the first dip is actually fixed by a single data point at ^ = 450 as can be seen 
from Fig. 3 in (de Bernardis et al. 2001). In order to remove these problems we approximate the 
experimental points in the range 250 < i < 850 (12 experimental points) by one single fifth order 
polynomial (6 free parameters). The number of degrees of freedom Nf = 6. The best-fit curve with 
xl^.^ = 3.26 is presented in Fig. 2. Its 6 coefficients are oq = 1.13 • 10^, ai = -944.7, 02 = 3.106, 
03 = -4.92 • 10-^, 04 = 3.75 • 10~^, 05 = -1.099 • lO"'^. This method allows us to determine the 
locations and amplitudes of both dips and of the second peak by taking into account the relatively 
prominent raises to the third and especially to the first peak. The local extrema of the polynomial 
best-fit give the following locations and amplitudes of the 1st and 2nd dip (positive curvature 
extrema) and the 2nd peak (negative curvature) between them. The la error bars are determined 
as above: 



4i 


= 4131^7, 




= 1960^272 






= 5441^^, 


A 


= 2266+27I 






= 746^^^^ 


M 


= I6O5+33I 





The results discussed here arc presented in Table 1 and shown in Figs. 1 and 2. In Fig. 2 also the 
la (Ax^ = 7.04), 2a (Ax^ = 12.8) and 3a (Ax^ = 20.1) confidence contours are shown. The 1-a 
contours for all peaks and dips are now closed. The 2a contour for the 2nd peak has a 'corridor' 
connecting it with the 3rd peak. As we have noted before (see Fig. 1), the 2a confidence contour for 
the 3rd peak is also open towards low i. This implies that we can not establish the second peak at 
2a C.L. The probability of its location is spread out over the entire range 450 < I < 920. Therefore 
at 2a C.L. we can not state whether the Boomerang-2001 results indicate a second peak without a 
third, or third without a second or both. We only can state at 2a C.L. that there are one or two 
negative curvature extrema of the function £{i+l)Ci/2Tr situated in the range 450 < i < 920 with 
amplitude in the range 1500 < £{£+l)Ci/2Tr < 2700/xif^. Furthermore, there is one secondary peak 
present in this i range at 2a confidence. Now, the contours for the dips are in logical agreement 
with the information about the second and third peaks. If at 2a C.L. the 2nd peak can be at the 
range of the location of 3rd one, then the 1st dip will move to £ ^ 520. Its 2a contour is closed 
since the 3rd peak has a closed 2a C.L. contour at the high £ side. On the contrary, the 2nd dip 
is open at high £ as it disappears when the 'second' peak disappears and the 'third' peak becomes 
the second. The 3a contours for the 2nd and 3rd peak as well as for 1st and 2nd dips are open in 
the direction of high £. 

So far we discussed only the statistical errors. The Boomerang LDB measurements have 
two systematic errors: 20% calibration uncertainty and beam width uncertainty leading to scale- 
dependent correlated uncertainties in the determination of the power spectrum (Netterfield et al. 
2001). The calibration error results in the same relative error for all data points and can be taken 
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into account easily. The beam width uncertainty which induces an error which becomes larger at 
higher values of ^, needs more care. 

We estimate the beam width uncertainty as follows: Using the data of Netterfield et al. (2001) 
for the la dispersion of the CMB power spectrum due to beam width uncertainty, we have to 
estimate its effect on determination of peak and dip locations and amplitudes. To take into account 
the effect of a la overestimated beam width we have lowered the central points of the CMB power 
spectrum presented in Table 3 of (Netterfield et al. 2001) by multiplying them by the ^-dependent 
factor 

f^(i) = 1 - 1.1326 • 10^"^^ - 2.72 • lO-'^^^ 

To take into account the effect of a la underestimated beam width we raise the central points of 
CMB power spectrum all by multiplication with the factor 

f^(i) = 1 - 6.99 • 10"^^ + 5.53 • 10-'^^^. 

(/o(^) and fu{^) are best fits to the Boomerang-2001 data and are shown in Fig. 3). For both cases 
we have repeated the peak and dip determination procedure. Best-fit values determined for the 
central points of CMB power spectrum give us the la errors of the peak/dip characteristics due 
to beam width uncertainties. The results are presented in Table 1. They show that error bars of 
all peak and dip locations caused by the beam width uncertainty are substantially less than the 
statistical errors. But they dominate for the amplitude of the 3rd peak and are comparable with 
the statistical error for the amplitudes of the 1st dip, the 2nd peak and the 2nd dip. However, the 
beam size errors are significantly smaller than statistical errors for the 1st acoustic peak. 

Since all sources of errors have different nature and are statistically independent they add in 
quadrature. The resulting symmetrized total errors are shown in the before last and last columns 
of Table 1. They are used in the cosmological parameter search procedure described in Section 4. 

2.2. Adding DASI and MAXIMA-1 data 

We have repeated the determination of peak and dip locations and amplitudes with the data of 
two other experiments, DASI (Halverson et al. 2001) and MAXIMA-1 (Lee et al. 2001), released 
simultaneously with Boomerang-2001. Both confirm the main features of the Boomerang CMB 
power spectrum: a dominant first acoustic peak at ^ ~ 200, DASI shows a second peak at ^ ~ 540 
and MAXIMA-1 exhibits mainly a 'third peak' at ^ ~ 840. The results presented in Figs. 4,5 and 
6 are quantitative figures of merit for their mutual agreement and/or disagreement. In Fig. 4 the 
1, 2 and 3a contours for the first peak location and amplitude for each experiment as well as the 
contours for the combined data are presented. 

In Fig. 4 one sees that MAXIMA, like Boomerang, indicates on the existence of the first acoustic 
peak at approximately 3a C.L. But its la contour for location of this peak in the + l)C^/27r) 

plane does not intersect the Boomerang la contour, though their projections on i and £{£+ 1)0^/21: 




Fig. 1. — Best-fit parabolas for the acoustic peaks of tlie Boomer ang-2001 CMB power spectrum 
{Xmin — 2-^9 fo'^ the 1st peak, Xmin — 0-92 '^^'^ peak and Xmin — 0-65 for the 3rd peak) 

as well as 1 (solid), 2 (dotted) and 3a (dashed) contours for their locations and amplitudes are 
shown. The crosses indicate the top of the best-fit parabolas. The contours limit the regions in the 
+ l)Ci/2TT) plane which contain the tops of parabolas with Ax^ = 3.53, 8.02, 14.2 for the 
1st and 3rd peaks (Nf = 3) and A^^ = 5.89, 11.3, 18.2 for the 2nd peak {Nf = 5). 
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Fig. 2. — The best polynomial fit for the Boomerang-2001 CMB power spectrum in the range of the 
1st dip, 2nd peak and 2nd dip {Xmin — 3-26), and the 1, 2 and 3a contours for their locations and 
amplitudes. The crosses indicate the positive (dips) and negative (peak) curvature extrema. The 
contours limit the regions in the + l)C^/27r) plane containing the corresponding extrema of 

polynomial fits with Ax^ = 7.04, 12.8, 20.1 {Nj = 6). 
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Fig. 3. — The correction factors /o(^) (upward pointing triangles) and /u(^) (downward pointing 
triangles) for the correlated CMB power spectrum error caused by the uncertainty of the effective 
beam width of the Boomerang experiment as given in Netterfield et al. (2001). The solid lines are 
the fitting functions /„ and /„ given in the text. 

axes do. This can be caused by systematic (normalization) errors inherent in both experiments. 

Approximately a quarter of the area outlined by the Boomerang 2a contour falls within the MAX- 
IMA 2(7 contour. The experiments show the same level of agreement in the data on 3rd acoustic 
peak (Fig. 5). In the range of the 1st dip - 2nd peak - 2nd dip, the MAXIMA data have no 
significant extrema, even la contours are open in both directions of the i axis. 

The DASI experiment establishes the location and amplitude of the first acoustic peak at 
somewhat more than la but less than 2a. The remarkable feature is the intersection of the la 
contours of DASI and Boomerang. Approximately 1/5 of the area outlined by the MAXIMA 2a 
contour is within the corresponding DASI contour. 

Our analysis has also shown that the DASI data on the second acoustic peak agree very well 
with Boomerang; the la contours nearly superimpose. The agreement of these two experiments is 
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Fig. 4. — The location of the first acoustic peak in the plane + l)C^/27r) for the Boomerang- 

2001 (blue diamond), DASI (green triangle), MAXIMA-1 (red square) data and for all experiments 
together (black cross) determined as maxima of corresponding best-fit parabola. The 1, 2, 3a 
confidence contours are also shown. Xmin^ Ax^ for the first peak of the Boomerang- 

2001 data are given in the caption of Fig. 1, for the other cases we have: DASI - Xmm ~ -'-■8' 
Nf = 1 and Ax^ = 1, 4, 9 for the 1, 2, 3a confidence contours accordingly, MAXIMA-1 - 
xl^^^ = 4.22, Nf = 2 and Ax^ = 2.3, 6.17, 11.8, all experiments together - xLn = 15.3, Nj = 9 
and Ax^ = 10.43, 17.18, 25.26. The dominant contribution to Xmin comes from the MAXIMA-1 
data. 



- 12 - 



< 



- XAII 

4000 -O Boomerang 
I □ MAXIMA- 1 ^ 




900 



950 



000 



Fig. 5. — The location of the third acoustic peak in the plane {£, i{l + l)Ce/2TT) for the Boomerang- 
2001 (blue diamond) and MAXIMA-1 (red square) data and for both experiments together (black 
cross), determined as maxima of the corresponding best-fit parabola. The 1, 2, 3cr confidence 
contours are also shown. Xmiw for the Boomerang-2001 third acoustic peak is given 



in the caption of Fig. 1, for the other cases we have: MAXIMA-1 - xi 



0.67, 



/ 



1 and 



Ax^ = 1, 4, 9, all experiments together - xl 



3.0, Nf = 6 and Ax^ = 7.04, 12.82, 20.06. The 



dominant contribution to Xmin comes from the MAXIMA-1 data. 
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impressive. 

We have repeated the determination of peak and dip locations and amplitudes using the data 
of all experiments jointly. The contours for the combined data are shown by the thicker black lines 
in Figs. 4 and 5. The C.L. contours for the 1st dip, 2nd peak and 2nd dip determined as regions 
of locations of negative and positive extrema of a 5-th order polynomial fit are shown in Fig. 6. 
The comparison with the corresponding figure for the Boomerang data alone (Fig. 2) shows their 
agreement. 

The best-fit values of ^p^, Ap^ {i = 1, 2, 3) and i^^,, A^^ {k = 1, 2) as well as their la statistical 
errors are given in Table 2. 

Other methods of model-independent determinations of acoustic oscillation extrema were pro- 
posed by (Douspis & Ferreira 2001; Miller & Nichol 2001). Their analysis like ours finds low sta- 
tistical significance (less than 2a) for the detection of second and third peaks. 

The peak locations and amplitudes from the Boomerang-2001 CMB data presented in the 
Table 1 show good quantitative agreement in the locations and, somewhat less good, in the am- 
plitudes obtained from the corresponding data of the other experiments and all the data together. 
The agreement can be improved when other error sources (calibration, beam width uncertainty, 
cosmic variance etc) of each experiment are taken into account. With some luck, the new mission 
MAP, which has been launched successfully last June, will remove many of the current problems 
and will considerably improve the data on CMB power spectrum. 

Clearly, the existence of the first peak in the spectrum is very well established in the present ex- 
periments. It has already been established (with less accuracy) before (see, eg. (Knox & Page 2000; 
Novosyadlyj et al. 2000)). Finally we have also assessed approximately the probability for each of 
the experiments to show no secondary peaks structure whatsoever. For the Boomerang data this 
probability is less than 4.6%, while for MAXIMA it is on the order of 15% and for DASI about 8%. 

3. Analytic determinations of the locations and amplitudes of the acoustic peaks 

and dips 

In order to use the data in Tables 1 and 2 to determine cosmological parameters, we need a fast 
algorithm to calculate the peak and dip positions for a given model. Here we improve the analytical 

approximations of peak/dip positions and amplitudes which have been derived in several papers 
(Efstathiou & Bond 1999; Hu et al. 2001; Durrer k Novosyadlyj 2001; Doran & Lilley 2001). We 
start by discussing the normalization procedure. 
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Table 1. Best fit values for locations {£p) and amplitudes {Ap, [/J^K"^]) of the peaks and dips in 
the CMB temperature fluctuation power spectrum measured by Boomerang (Netterfield et al. 
2001). Statistical errors (1st upper/lower values) and errors caused by beam width uncertainties 
(2nd upper /lower values) are shown in columns 2 and 3. The 20% calibration uncertainty is 
included in the symmetrized total errors presented in the last column. 



Features ip Ap Ip Ap 



1st peak 


919+13+2 
^^^-20-3 


542fi+540+112 
>J^^O_539_136 


212 ± 17 


5426 ± 1218 


1st dip 


4,0+54+6 


iqRn+2''2+142 
iyDU_282-158 


413 ± 50 


1960 ± 503 


2nd peak 


544+56+14 
0*^-52-14 


99fiR+275+309 
zzou_274_283 


544 ± 56 


2266 ± 607 


2nd dip 


746+114+9 


lfin^+3'^3+422 


746 ± 89 


1605 ± 650 


3rd peak 


8431^^15 


9077+426+720 
' _4li_573 


843 ± 35 


2077 ± 876 



Table 2. Best fit values for the locations and amplitudes of peaks and dips in the CMB 
temperature fluctuation power spectrum from the DASI and MAXIMA-1 experiments. Statistical 
errors are determined as described in the text. In the last column the results obtained from the 
data of all three experiments together are presented. 



DASI MAXIMA-l All three experiments 



Features 
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y4.p 




j4.p 




Ap 


1st peak 
1st dip 
2nd peak ^ 
2nd dip 
3rd peak 


378115 
5361^0 

709l« 


471611^5 

157811?" 
236211?^ 
1799l|gi 


2361^° 

475llf 
435 - 739'' 
435 - 739 
813l??l 


4438l?« 

1596lS^ 
1500 - 2800 ^ 
1000 - 2700 ^ 

^o2»_-^5g4 


40611^ 
545ir 
736115? 

»47+252 
''*'-146 


5041+10" 

18431385 
22661^9^ 
16611^1^ 
2175l«^^ 



''The cxtrcma were determined by approximating the experimental CMB power spectrum 
by 5-th order polynomial 

''Just the ranges where the probability to find the peak or dip is > 68.3% are indicated 
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Fig. 6. — The locations of the second acoustic peak and the two dips in the {£,£{£ + l)Ci/2Tr) 
plane are shown together with the best-fit 5-th order polynomial. All points in this range from 
Boomerang, MAXIMA-1 and DASI have been used jointly. The second peak and the dips are 
determined as extrema of negative and positive curvature of the corresponding best fit polynomial. 
The 1, 2, 3a confidence contours are also shown. For all experiments together, the best-fit gives 
Xmin — ^'^■^ — 1^ degrees of freedom. 
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3.1. Normalization of the density power spectrum 

The 4-year COBE data, which estabhsh the amphtude and the form of the CMB power spec- 
trum at the largest angular scales (i < 20), are taken into account via the approximation for Cio 
proposed by Bunn & White (1997). This requires accurate calculations of Q in the range i < 12. 
The dominant contribution on these angular scales is given by ordinary Sachs- Wolfe (SW) effect. 
However, the Doppler (D) effect and the cross-correlation term Sachs- Wolfe - adiabatic (SW-A) in 
the general expression for the correlation function (^(ni) • ^(n2)) have to be taken into account 
as well if we want to achieve an accuracy better than 20% (see Appendix A). For A dark matter 
models and models with non-zero 3-curvature, also the integrated Sachs- Wolfe effect (ISW) con- 
tributes. We use the factors (> 1) introduced and calculated by Kofman & Starobinsky (1985) 
and improved by Apunevych & Novosyadlyj (2000), so that cf^^^^^ = KjCf^ (for details see 
Appendix A). 

The normalization of the power spectrum of scalar perturbations then consist in two steps: 

i) We calculate 

Q = CSW+ISW ^ ^ ^ ^SW-A 

(for i = 2, 3, 5, 7, 10, 11, 12) by the analytical formulae given in Appendix A with arbitrary 
normalization. This determines the shape of the CMB power spectrum in the range of the COBE 
data, and hence the best-fit parameter C^^^^ to 4-year COBE and the first and second derivatives 
as defined in Bunn & White (1997) for models with given cosmological parameters; 

ii) Since each term in the expression (1) is oc (5^, where 6h is the present matter density pertur- 
bations at horizon scale, we can now determine and along with it the value of the normalization 
constant for scalar perturbations Ag = 27r^(5^(3000Mpc//i)^+"« for a model with given cosmological 
parameters. Here Ug is the spectral index for primordial scalar density perturbations and h is 
dimensionless Hubble parameter (in units of 100 km/sec/Mpc). 

Both these steps are also performed in CMBfast. Hence, our normalization procedure for the 
power spectrum is equivalent to normalization with CMBfast. Calculations show that our value 
C}q-*^^ never differs from the result of CMBfast by more than 3%. The accuracy of the overall 
normalization constant 6^ for ADM models with appropriate values of parameters is better then 
5%. This has been controlled by comparing the value of as from CMBfast with our semi-analytical 
approach. This error simply reflects the accuracy of the analytical approximation of the transfer 
function for density fluctuations by Eisenstein &; Hu (1999) which we have used. 

3.2. Positions and amplitudes of CMB extrema: analytic approach 

One of the main ingredients for our search procedure is a fast and accurate calculation of the 
positions and amplitudes of the acoustic peaks and dips, which depend on cosmological parameters. 
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The dependence of the position and amplitude of the first acoustic peak of the CMB power 
spectrum on cosmological parameters has been investigated using CMBfast. As expected, the 
results are, within reasonable accuracy, independent of the hot dark matter contribution (Oj/). 
This was also shown by Novosyadlyj et al. (2000). For the remaining parameters, Us, h, fib, ficdm 
and fl\, we determine the resulting values ip-^ and Ap-^ using the analytical approximation given by 
Efstathiou & Bond (1999) and Durrer &; Novosyadlyj (2001). In these papers the CMB anisotropy 
spectrum is approximated in the vicinity of the first acoustic peak by 



27r 27r ^ ^ ^ 



A{^b, ^cdm, ^k, ris, h) exp 



(2) 



where Aip^ = 0A2£p^ , Cf^ is the Sachs- Wolfe approximation for the Qs derived in the Appendix 
A, Eq. (A9), and 

exp [ai + a2UJcdm + a^^cdm + aAUJh + 05^6+ 

+aGU}b(^cdm + aruJk + agujl + ^ding - 1)]- (3) 

Here cob = f^bh'^, ^cdm = flcdmh"^, = (1 — f^m — ^^a)^^- The position of the acoustic peaks is deter- 
mined as in (Efstathiou &: Bond 1999) for open and flat models and in (Durrer &: Novosyadlyj 2001) 
for closed models. The coefficients are defined by fitting to the numerical CMBfast amplitudes 
of the first acoustic peak on a sufficiently wide grid of parameters. We find: ai = 2.503, 02 = 8.906, 
03 = -7.733, 04 = -115.6, 05 = 35.66, = -7.225, 07 = 1.96, ag = -11.16, ag = 4.439. The accu- 
racy of the approximation is better than 5% in the parameter range 0.2 < < 1-2, < ft\ < 0.8, 
0.015 < < 0.12, 0.8 < < 1.2 and 0.4 < h < 1.0. The approximation for the amplitude 
breaks down in the models with large curvature (O^ < —0.2 and fl^ > 0.6) and low baryon density, 
{ub <. 0.006). 

To calculate the amplitudes of the 2nd and 3rd peaks, we use the analytic relations for the 
relative heights of these peaks w.r.t the first peak as given by Hu et al. (2001). 

= Ap-iH2{flmjflbi''T's)j 

Ap3 = ^PiH3{^m,^b,ns), (4) 

where the functions H2 and H3 are given in Eqs. (B16) and (B17) respectively. For the locations 
of 2nd and 3rd peaks we use the analytic approximations given by Hu et al. (2001) and Doran &i 
Lilley (2001) (see also Appendix B). 

Unfortunately, we have no analytic approximation for the dip amplitudes and hence we can 
not use their experimental values to determine cosmological parameters. But a sufficiently accurate 
analytic approximation for the location of the 1st dip is given in Doran &; Lilley (2001). We use 
it here. 
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Hence, we have analytical approximations for the dependences of the positions and amplitudes 
of three acoustic peaks and the location of the 1st dip on cosmological parameters 

(■Pii^m, ^A, ^k, ^b, h), 

Ap^{nm,^A,^k,^b,ns,h), {i = 1,2,3) 

^di {^m, ^A, ^6) h). 

Comparing the analytical values for different sets of parameters with numerical calculations using 
CMBfast, shows that the accuracy is about 5% for all locations and amplitudes of 1st and 3rd 
peaks in the ranges of cosmological parameters indicated above. The accuracy for amplitude of the 
2nd peak is always better than 9% in the same ranges. For some parameter values the second peak 
is underestimated. 

Of course, in principle the half- widths of peaks and dips, their "smoothness", "sharpness" or 
their inflection points may contain additional information on cosmological parameters. But when 
one compares CMB spectra obtained in the different experiments up to date, one concludes that 
their accuracy is not sufficient to influence the resulting cosmological parameters at the present 
precision. Therefore, we only use the most prominent observable patterns of the CMB power 
spectrum - locations and amplitudes of acoustic peaks. 

For the convenience, we present all analytic approximations used here in Appendix B. 



4. Cosmological parameters from the CMB peak amplitudes and locations 

We now use the results of Sections 2 and 3 to determine the cosmological parameters Qrn, ^A, 
Oj, (one sort of massive neutrino), Q^, Us, h, T/S (= d^"-^"^ / CIq"''''^^) and Tc (optical depth to de- 
coupling). We use the method described in detail in a previous paper (Durrer & Novosyadlyj 2001). 
We include 8 experimental points from the CMB power spectrum (COBE Cio, the amplitudes and 
locations of three acoustic peaks and the location of the first dip). In order to have a positive num- 
ber of degrees of freedom, Nf > 1 we add a weak constraint for the Hubble constant, 0.5 < h < 0.8 
when searching for 8 parameters. 

In order to establish la confidence intervals for each parameter we have applied the marginal- 
ization procedure described in (Durrer & Novosyadlyj 2001). The results are presented in Table 3. 

At first we check how the best fit parameters depend on the accuracy of the peak/dip locations 
and amplitudes. For this we compare the resulting cosmological parameters from the data given in 
Table 1 if we consider only statistical errors and with total errors. The results are given in the first 
two rows of Table 3. In spite of the different errors of the experimental values and their relations 
(the total errors of peak amplitudes increase faster with the peak number than the statistical error) 
the best fit cosmological parameters are similar. 

In order to estimate the sensitivity of cosmological parameters to experimental values we 
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substitute the Boomerang data on peak/dip locations and amplitudes by the values from Table 2, 
obtained from all three experiments combined. The comparison of the results in the third row with 
those above shows that the results are practically unchanged. We believe that the Boomerang data 
on peak/dip locations and amplitudes are best studied and have well established statistics, hence 
we use them in the following determinations. 

The neutrino contents for these three data sets can be rather large, ~ 0.3. This is due 
to the low sensitivity of the CMB anisotropy spectrum to Q^. When il,^ = is fixed, the best-fit 
values for the remaining parameters stay practically unchanged and also increases only very 
little. The CMB doesn't care whether dark matter should be hot or cold. In order to distinguish 
between cold and hot dark matter data which is sensitive to the density power spectrum on smaller 
scales needs to be added. In Table 3 we therefore exclude from the determination procedure 
and fix its value to in the first five rows. Contrary, for T/S and Tc we obtain 0, but with large 
la confidence limits due to the degeneracy in T/S, Tc and (Efstathiou & Bond 1999). 

A remarkable result is the good agreement of the best-fit content of baryons ^i,h^ ~ 0.02 
with the constraint from standard nucleosynthesis and the observed intergalactic content of light 
elements (Buries et al. 2001). The large la confidence limits for 0,^ and are due to the well 
known degeneracy of the CMB power spectrum in these parameters (Efstathiou & Bond 1999). 
However, the sum of their best fit values, $7^ -|- JIa = 1 — $7^ is always very close to 1 which implies 
that spatial curvature is small for the best fit model. 

We repeat the search procedure for different combinations of the CMB power spectrum extrema 
data (Table 1) with other cosmological data sets. The LSS data set used here ranges from the Lyman 
alpha forest, determining amplitude and spectral index of the matter power spectrum at very small 
scales, to large scale bulk velocities, cluster abundances and Abell cluster catalogs which determine 
as and the position of the 'knee' in the matter power spectrum. All of this is extensively discussed 
in (Durrer & Novosyadlyj 2001). 

The results for cosmological parameters from different combinations of observational data are 
shown in the lines 4 to 8 of Table 3. 

Adding a stronger constraint on the Hubble parameter, h = 0.65 it 0.10, and the big bang 
nucleosynthesis (BBN) constraint changes the best fit cosmological parameters only slightly. The 
SNIa constraint on the relation between JIa and (Perlmutter et al. 1999) substantially reduces 
the errors of these parameters (5th line in Table 3) as it removes the degeneracy between them. 
This degeneracy is also removed when we combine CMB and LSS data. 

The cosmological parameters obtained from the Boomerang CMB power spectrum extrema 
data combined with all other cosmological measurements (a detailed list can be found in Durrer & 
Novosyadlyj 2001) are presented in lines 6th, 7th and 8th of Table 3. The best fit values for the 
tensor mode amplitude T/S defined as CJq/C'iq in the last two cases are practically zero but the la 
confidence limits are wide due to the degeneracy of the CMB extrema in Ug, T/S and Tc- Even when 
combining the CMB with LSS data, the degeneracy in and T/S is not significantly removed. 
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Table 3. Cosmological parameters from the extrema of the CMB angular power spectrum in 
combination with other cosmological data sets. The upper/low values show la confidence limits 
which are obtained by maximizing the (Gaussian) 68 percent confidence contours over all other 
parameters. The LSS data set is the same as in (Durrer &: Novosyadlyj 2001). 
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*)This pairameter is fixed to 0. 

*The big bang nucleosynthesis constraint on baryon content fJjj/i^ = 0.02 ± 0.001 from Buries et al. (2001) is included. 

''The constraint on the Ha — relation from SNIa distance measurements (Perlmutter et al. 1999), [Clm — 0.75f2A] = —0.25 ± 0.125 
is added. 

'^In addition to the parameters given in the different columns, ■we have also to determine the Abell-ACO biasing parameter, bd- The 
result is: = 2.64 ± 0.27 

"^For this data set we obtain bd = 2.47 ± 0.19 

"bci =2.5 ±0.2 
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This is so, since a blue spectrum, which allows for a high tensor contribution to the CMB, can be 
compensated with a neutrino component which leads to damping of the matter power spectrum on 
small scales. If massive neutrinos are not allowed, the degeneracy between Ug and T/S is lifted as 
soon as small scale LSS data is included. 

The best-fit values of spectral index Ug in all cases are in the la range of the value obtained 
from the COBE 4-ycar data, = 1.2 ± 0.3 (Bennett et al. 1996; Gorski et al. 1996). When using 
the best fit model to calculate the data used to find it, practically all results are within the la error 
range of the corresponding experimental data. Only two out of 31 experimental points are slightly 
outside. Namely the best-fit value of — 0-75O.\ in the last determination is at l.la lower of its 
experimental value followed from SNIa test and as constraint established by (Bahcall Sz Fan 1998) 
from the existence of three massive clusters of galaxies is at l.lScr higher than model predicted 
value. But the value of as in our best-fit model, as = 0.91, is in the range of current estimates 
from the Sloan Digital Sky Survey a^^^^'^^^^ = 0.915 ± 0.06, (SDSS Collaboration 2001)), which 
is not included in our data set. The high degree of consistency within completely independent 
cosmological data sets is very encouraging. 

Moreover, all parameters of our best-fit model agree well with those extracted from the full 
Boomerang data (Netterfield et al. 2001) combined with LSS and SNIa priors (compare our values 
derived from CMB-I-LSS in the 7th row of Table 3 and theirs in the 4th row of Table 4). But our 
la ranges for most parameters are significantly wider. That is due to two facts. First of all, we 
allow also for a tensor component and neutrinos. This increase in the number of parameters also 
increases the degeneracies (e.g. between the tensor amplitude and the spectral index n^) thereby 
enlarging the errors of the physical parameters. In this sense our parameter errors arc rather 
to be compared with those of Wang et al. (2001) which allow for roughly the same degrees of 
freedom. But even their parameter estimation is somewhat more precise than ours. This is because 
we use more conservative errors for the peak and dip locations and amplitudes which include 
statistical, normalization and beam uncertainties. We think that with present data, and with the 
model assumptions made by our choice of parameters, our precision is realistic. Clearly, future 
experiments like the MAP satellite will improve this situation. 

Finally, to compare with the cosmological parameters obtained in our previous paper (Durrer 
& Novosyadlyj 2001, Table 4), where the same LSS data set was used, we have repeated the search 
procedure fixing T/S = Tc = 0. The best-fit values of the parameters with la errors obtained by 
maximizing the confidence contours over all other parameters are given in the last row of Table 
3. Comparing them with values in the last column of Table 4 from (Durrer Sz Novosyadlyj 2001) 
shows that both determinations have best-fit values in the la confidence limits of each other. There 
the best fit model has a slight positive curvature, here a slightly negative. The la confidence ranges 
here are somewhat wider than those obtained in the previous determination. These differences are 
due to the different CMB observable data set and the different normalization procedure. Even 
though in our previous analysis we have only taken into account the first peak. The errors in its 
location and amplitude were significantly underestimated, leading to smaller error bars. 
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5. Conclusions 

We have carried out a model-independent analysis of recent CMB power spectrum measure- 
ments in the Boomerang (Netterfield et al. 2001), DASI (Halvcrson et al. 2001) and MAXIMA 
(Lee et al. 2001) experiments and we have determined the locations and amplitudes of the first 
three acoustic peaks and two dips as well as their confidential levels (Table 1-2, Fig. 1-7). 

In the Boomerang experiment the second and third acoustic peaks are determined at a confi- 
dence level somewhat higher than la. Experimental errors which include statistics and systematics 
are still too large to establish the secondary peak locations and amplitudes at 2a C.L. Only the 
position of one (the third ) secondary peak can be bounded from above < 900, at 2a C.L. The 
same situation is encountered when determining the locations and amplitudes of the first and second 
dips. However, the location and amplitude of the first peak, are well established with confidence 
level, higher than 3(T. 

The MAXIMA experiment also shows the existence of the first acoustic peak at approximately 
the same confidence level as Boomerang. But the \a contours for the peak position in the plane 
{£,£(£ + 1)C£/2tt) do not intersect. However, their projections onto the £ and £(i + 1)C£/2tt axes 
do. Approximately one quarter of the area inside the Boomerang 2a contour falls within the 
corresponding MAXIMA contour. The same level of agreement of these experiments is found in 
the data on the 3rd acoustic peak. In the range of 1st dip - 2nd peak - 2nd dip the MAXIMA data 
give no significant information. Even the la contours are open in both directions of the £ axis. 

The DASI experiment establishes the location and amplitude of the first acoustic peak at 
somewhat higher than la C.L. but less than 2a. The la contours for the position of the first peak 
of the DASI and Boomerang experiments intersect. Approximately 1/5 of the area outlined by the 
MAXIMA 2a contour is within the DASI 2a contour. The DASI data on second acoustic peak is 
in excellent agreement with the Boomerang results - the la contours practically coincide. 

We have also determined the locations and amplitudes of the acoustic peaks and dips using 
the data of all three experiments. The results are very close to those from the Boomerang data 
alone. 

To determine cosmological parameters from these data, we have improved the analytical ap- 
proximations for the peak positions and amplitudes to an accuracy (determined by comparing the 
approximations with the results of CMBfast) better than 5% in a sufficiently wide range of parame- 
ters. Wc have also developed a fast and accurate analytical method to normalize the power spectrum 
to the 4-year COBE data on Ciq. Our analytical approximation is accurate to a few percent (in 
comparison to CMBfast) when all main effects (ordinary Sachs- Wolfe effect, integrated Sachs- Wolfe 
effect, adiabatic term, Doppler term and their mutual cross-correlations) are taken into account. 
For example, in the model with parameters presented in the last row of Table 3 the relation of con- 
tribution from these components at ^ = 10 are Cf^ : : Cf^~^ : Cg, = 1 : 0.098 : -0.24 : 0.42. 

The cosmological parameters extracted from the data on locations and amplitudes of the first 
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three peaks and the location of the first dip are in good agreement with other determinations 
(Netterfield et al. 2001; de Bernardis et al. 2001; Pryke et al. 2001; Wang et al. 2001; Durrer & 
Novosyadlyj 2001). That shows also that present CMB data can essentially be compressed into the 
height and slope of the Sachs- Wolfe plateau (at i = 10) and the positions and amplitudes of the 
first three acoustic peaks and the first two dips. 

A remarkable feature is the coincidence of the baryon content obtained from the CMB data, 
ilfe/i^ « 0.02 with the value from standard nucleosynthesis (0.02 it 0.001) (Buries et al. 2001). 
Moreover, the CMB data together with constraints from direct measurements of the Hubble con- 
stant, the SNIa data, the baryon content and the large scale structure of the Universe (the power 
spectrum of rich clusters, the cluster mass function, the peculiar velocity field of galaxies, Lj-a ab- 
sorption lines as seen in quasar spectra) select a best-fit model which gives predictions within about 
la error bars of all measurements. The cosmological parameters of this model arc Q\ = 0.64:~^q'2^, 
Qrn = 0.36l{]:^i, Qb = 0.047i[j:{]^^, Us = l-Otoir^ ^ = O-SSIoJt and = 0.15l{]:?^. The best-fit 
values of and T/S are close to zero, their la upper limits are < 0.17, T/S < 1.7. 

The cosmological parameters determined from the CMB acoustic peak/dip locations and ampli- 
tudes data show good agreement with other cosmological measurements and indicate the existence 
of a simple (adiabatic) best-fit model for all the discussed cosmological data within the accuracy 
of present experiments. 

It is a pleasure to acknowledge stimulating discussions with Alessandro Melchiorri, Roman 
Juszkiewicz and Roberto Trotta. BN is grateful to the Tomalla foundation for a visiting grant and 
to Geneva University for hospitality. RD thanks the Institute for Advanced Study for hospitality 
and acknowledges support from the Monell Foundation. 

APPENDIX 



A. An analytic approximation for the CMB power spectrum at large angular scales 

On sufficiently large angular scales (larger than the Silk damping scale) temperature fluctua- 
tions in the CMB can be related to density, velocity and metric perturbations at the last scattering 
surface and at later times by integrating the geodesic equation, similar to the classical paper by 

Sachs and Wolfe (1967). Here we discuss only scalar perturbations. Tensor perturbations can be 
simply added to the result and do not pose any significant difficulty. Scalar perturbations gener- 
ate CMB temperature fluctuations which can be written in gauge-invariant form as a sum of four 
terms - the ordinary Sachs- Wolfe effect, the integrated Sachs- Wolfe term, the Doppler term and 
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the acoustic term (Durrer 1990). 

DriVdec^dec) + ^i(%ec, Xdec)"-* + - 'i>)(,Vdec,^dec) 

rm 

■/ -^'){r^,^{r^))dri . (Al) 

Here $ and ^' are the Bardeen potentials (Bardeen 1980), Vi is the baryon velocity and is a gauge 
invariant variable for the radiation density fluctuations. A prime denotes the partial derivative 
w.r.t. conformal time 77. For perfect fluids and for dust we have ^ = — In Newtonian limit the 
Bardeen potentials just reduce to the ordinary Newtonian potential. For adiabatic perturbations 
\Dt. = \8„i — 1$ (see e.g. Durrer & Straumann (1999)), where 8,n is the usual matter density 
perturbation., it corresponds to in Bardeen's notation. The variables rj and x are conformal 
time and comoving position. 

In realistic models, cosmological recombination and decoupling of radiation from matter take 
place when pm > Pr- Hence the large angular scale CMB power spectrum can be expressed in 
the terms of solutions of Einstein's equations for adiabatic linear perturbations in a dust Universe. 
The CMB anisotropies on angular scales > 10" {£ < 20) are generated mainly by the linear 
perturbations of matter density, velocity and the gravitational potential at scales much larger than 
the particle horizon at decoupling. Our approximation makes use of these facts. 

We use the solutions of Einstein's equations for linear density perturbations in flat models of a 
Universe with dust and a cosmological constant which can be found in (Kofman &; Starobinsky 1985; 
Apunevych & Novosyadlyj 2000). The growing mode of density, velocity and gravitational potential 
perturbations, using the gauge- invariant variables introduced by Bardeen (1980) and normalizing 
the scale factor a{to) = oq = 1, are 

C{k) is (up to the time dependent factor Ks) the Fourier transform of the Bardeen potential, 
so that $(t,x) = (27r)-3/2 J $(t, k)e*i'''d3/c, 5™(t,x) = (27r)-3/2 J 5{t,\i)e}^'^d^k and F"(t,x) = 
(27r)-3/2 / \L)e'^^d^k. The factors Ks{t) = | (l - d/a^ adt^ and Kv{t) = | (d/a^ - d/ad) /g adt 
are both in the range < K, < 1 and reflect the reduction of growth of perturbations caused by 
the cosmological constant. The scale factor of the background model is given by 

= [t^J ^^"'^ [ — 2 — ) 

Here Hq = (d/a)(to) = d(to) is the Hubble constant today. The ii'.-factors go to 1 when t <^ to 
or when 0,^ 1 (^A ^0). At decoupling Ks = Ky = 1. An analytical approximation for 
Ks{t)/Q.rn with sufficient accuracy can be found in (Carroll et al. 1992) and for Ky(io)/^m in 
(Lahav et al. 1991). 



( ^ I (%,xo,n) = 
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The power spectrum of density fluctuations is given by 

P{k,t) = <5{t,k)6*it,k)>=Ask''^Tl{k;t)a\t)Klit)/nl, (A3) 
As = 27r2j^(3000Mpc//i)3+'*% 

where Tjn{k,t) is transfer function (divided by the growth factor) and Sh is the present matter 
density perturbation at horizon scale. We use the analytical approximation of Tm{k,t) in the space 
of cosmological parameters h, fi^, fib, ^A, and N^, (number of species of massive neutrino) by 
Eisenstein & Hu (1999). 

Prom Eq. (Al), taking into account adiabaticity and setting xq = 0, we obtain 

^(n) = l^iVdec, n??o) + 2 / ^^('^o ^'^^ du + naV{r]dec, n%) + l^mivdec: nr?o), (A4) 
T 3 Jo dr) 3 

where n is the unit vector in direction of the incoming photon and we have used x(r7) = 11(770 — ''7) , 
Xdec — nr/o- The variable u is the affinc parameter along the geodesic which begins at the observer 
and ends in the emission point at the last scattering surface. The present value of conformal times, 
rjQ gives also the present particle horizon or the distance to the last-scattering surface. The first 
term in (A4) is the well known Sachs- Wolfe effect (SW), the second term is the integrated Sachs- 
Wolfe effect (ISW) which is important only at late times, where Ks{t) starts to deviate from 1 and 
|i| 7^ 0, the third is the Doppler term (D) and the last is the acoustic term (A). At large angular 
scales (~ 10"), where anisotropics have been measured by COBE (Bennett et al. 1996), the SW 
and ISW effects dominate. However, if we want to calculate Cio with good accuracy, we must to 
also take into account the other terms. The angular correlation function of AT/T can be written 
symbolically as 

AT AT 

< — (ni) • — (na) >= < SW • SW > +2 < SW • ISW > + < ISW • ISW > 

+ <AA>+2<SWA> + <DD>. (A5) 

The cross-correlators < D • SW > and < D ■ A > are omitted because they are strongly suppressed 
on large angular scales. Indeed, if one uses Fourier presentations for the variables (A2) in the 
equations (A4-A5) one finds that the A;-integrand of these terms contains a spherical Bessel function 
ji(kriQ(ni — n2)) which oscillates for large angular separations, strongly reducing the integral if 
compared to the < SW • A > term where the integrand has a definite sign. The terms < ISW • A > 
and < ISW • D > are also omitted because the ISW effect gives the maximal contribution to AT/T 
at the largest angular scales of the range of interest (at lowest spherical harmonics) where A and 
D are nearly zero. At 'smaller' angular scales {£ ^ 10) where contribution of A and D are not 
negligible, the ISW effect is very small. Therefore, their cross-correlation terms are very small. 

We develop the n-dependence of ^ (n) in spherical harmonics 

AT s-^ 

^('^) = 2^a«m(f?o)y^m(n), < aemal'm' >= SemSt'm'Ce . 
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The CMB power spectrum, C^, has the same components as the correlation function: 

Q = cf w + c|w-isw + ^ISW + + ^SW-A ^ (^6) 

Each component on the right hand side comes from the corresponding contribution to AT above 
and is proportional to 5^. Using the solutions (A2) we obtain analytic approximations for them. 

We first approximate the SW and ISW contributions in the form 

^SW+ISW _ ^SW _|_ ^SW-ISW _|_ ^ISW ^ j^2qSW 

where the factors Kg^ (> 1) take into account the contribution of the ISW effect for each spher- 
ical harmonic. They have been calculated by Kofman k. Starobinsky (1985) and Apunevych & 
Novosyadlyi (2000) for different A-models. Instead of the direct time consuming calculations of the 
ISW contribution, we use the following analytic approximations: 

Kl = l + 8.20423 X exp(-fi^/0.01157) + 3.75518 x exp(-fi^/0.13073), 
Kl = l + 2.25571 X exp(-J^^/0.03115) + 2.35403 x exp(-J^„/0.15805), 
Kl = l + 1.80309 X exp(-J^^/0.0323) + 1.88325 x exp(-n^/0.16163) 
and Kj = l + [23.46523 x exp(-J7^/0.0122) + 11.03227 x exp(-J^^/0.14558)]/(^ + 0.5) 

for i > 5. These approximation formulae are determined from the data presented in the tables of 
(Kofman & Starobinsky 1985) and (Apunevych & Novosyadlyj 2000). 

Using solutions (A2) and the definition of the density power spectrum (A4) we obtain the 
following general expression for the SW contribution to the CMB power spectrum: 

= 2n!-'w^{to) i dkk^'~'Tl{t,e,, k)jnkm), (A8) 

where T^{tdec,k) is the transfer function of matter density perturbations at decoupling, D{to) = 
Ks{to)/flm is the value of the growth factor at the current epoch, and is the spherical Bessel 
function of order £. For reasonable values of spectral index — 3 < < 3 the main contribution to 
the integral (A8) comes from very small k where T^itdec: k) 1 and can be omitted. Then integral 
can be performed analytically and the result can be expressed in terms of F-functions: 

^sw ^ ^'^1 r(3-n.)F(£ + V) 

^ 8D^{to)r'^{2-ns/2)r{e + ^y ^ ' 

In the same way we obtain the expressions for the other components of equation (A6): 

CfW-^ = - /o° i ml to / dkk-'T^itdec, k)n{tdec, fc)j|(fe%), (All) 

_„ns+lr2 ('j. \ POO 

^? = pyiS;i <iM»-T.^(W, (A12) 
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where Ti,{t(iec) is transfer function for density perturbations of baryons (Eisenstein & Hu 1998) and 
(') is the derivative w.r.t the argument x = kr)o. The minus sign in the expression for Cf^~^ 
reflects the anti-correlation of the gravitational potential and density fluctuations: large positive 
density fluctuations generate deep negative potential wells. If wc set = Tf, = 1, the integrals 
(All -A12) diverge for all i for > 1 because the main contribution to integrals of the D and A 
terms comes from small scales. Hence here the transfer functions must be kept and the integrals 
have to be calculated numerically. Fortunately, the integrands decay rapidly for large wave numbers 
and 99.9% of the contribution comes from the range 0.001 < A;%ec < 0.1, so that the integration is 
not very time consuming. 

In Fig. 7 the CMB power spectrum Q at large angular scales {9 > 20", £ < 20) together with 
the contributions from the different terms given in (A9-A12) is shown for a pure matter and a A— 
dominated model. The relation of the contributions from different terms at ^ = 10 are 

CfW : : Cf^'-^ -Cf = 1: 0.04 : -0.11 : 0.22 

for the matter dominated flat model (17^ = 1) and 

CfW : : Cf^-^ : = 1 : 0.08 : -0.23 : 0.39. 

for the A dominated model with the cosmological parameters shown in the figure. Therefore, a few 
percent accuracy of the normalization to 4-year COBE Cio data can be achieved only if all these 
effects are taken into account. 

In Fig. 8 the CMB power spectrum at large scales calculated using the analytic formulae 
(A9-A12) and using CMBfast are shown for comparison. In the left panel we also present the 
power spectrum calculated by the analytical approach of (Hu &; Sugiyama 1995) (renormalized to 
the CMBfast value of Cio). 

The calculations show that value of Cio calculated by our method deviates from the value 
calculated with CMBfast by 0.5% for the matter dominated flat model (17^ = 1) and 2.7% for the 
A dominated model (O^ = 0.2). Therefore, our analytic approach is sufficient to normalize fast the 
power spectrum of scalar perturbations to the 4-year COBE data with virtually the same precision 
as CMBfast, the difference is less than 3%. (Remember, that the experimental errors of the COBE 
data are about 14%, so that the best-fit normalization parameter Ci(f^^ has the same error.) 

Another comparison of our normalization procedure with CMBfast comes from the value of as- 
For the fiat model (left panel of Fig. 8) our approximation for the normalization together with the 

analytical transfer function of Eisenstein & Hu (1998; 1999) leads to cjg = 1.58, the corresponding 
value calculated from CMBfast is 1.53. For the A dark matter model (right panel) our erg = 0.62, 
while CMBfast gives as = 0.64. The agreement of both approaches is quite well (the 5% difference 
includes also the errors in the approximation of the transfer function which is actually of this order) . 

The slight deviation of Cio as calculated by our code from the value obtained with CMBfast 
(< 3%) in spite of using the same analytic best-fit formula for C^^^^^ by (Bunn & White 1997) 
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Fig. 7. — The CMB power spectrum and the different contributions discussed in the text (formulae 
A9-A12) for a pure matter model (left panel) and a A dominated model (right panel). 




Fig. 8. - The CMB power spectrum at COBE scales calculated by CMBfast (solid line) and by our 
analytical formulae (A9-A12) (dotted line). For the pure matter model we also show the spectrum 
calculated with the analytic approach of (Hu & Sugiyama 1995) (dashed line) (this approach does 
not allow a cosmological constant). All spectra are normalized to the best fit for Cio from the 
4-year COBE data given in (Bunn & White 1997). 
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is due to a difference in the form of the spectra as shown in the Fig. 8. This difference grows 
when rijji decreases. There are several possible reasons for this deviation in the form of the CMB 
power spectrum in our analytic approach from the exact numerical calculation: 1) We have used 
the solutions for the evolution of density, velocity and gravitational potential perturbations in 
the A- dust Universe. In reality, at decoupling the role of radiation is not completely negligible, 
this slightly influences the dynamics of the scale factor and the evolution of perturbations. It also 
results in additional time dependence of gravitational potential (early integrated Sachs- Wolfe effect) 
which is not taken into account here. 2) Our approach does not take into account the effects of the 
collisionless dynamics of photons and neutrinos after decoupling. Especially, the induced anisotropic 
stresses lead to ~ 10% difference of the gravitational potentials in the radiation-dominated epoch 
which results into a corrections of a few percent in the C^'s. 3) Instantaneous recombination 
and tight coupling which were assumed, also cause slight inaccuracies. They should, however, be 
extremely small on the angular scales considered here. 4) To calculate the terms Cf, Cf'^"^ and 
we have used the analytic approximations for the transfer functions, Tjn{tdec, k) and Tb{tdec, k), 
by Eisenstein &; Hu (1998; 1999) which have an accuracy ~ 5%. 

More details on the theory of CMB anisotropics can be found in the reviews by Durrer & 
Straumann (1999) and Durrer (2001). 
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B. Analytic formulae for the amplitudes and locations of acoustic peaks and dips in 

the CMB power spectrum 

For completeness, we repeat here the formulas used in our parameter search which can also be 
found in the cited literature. 

We assume the standard recombination history and define the redshift of decoupling z^i^^ 
the redshift at which the optical depth of Thompson scattering is unity. A useful fitting formula 
for Zdec is given by (Hu &: Sugiyama 1996): 

Zdec = 1048[1 + 0.00124u;-0-^38]^;^ ^ ^^^^2]^ (Bl) 

where 

gi = 0.0783a;j-°-^^^[l + 39.5a;^-^^^]-\ 52 = 0.56[1 + 21.1ujI-^'^]-^ , 
LOb = Ofe/i^ and ujm = ^mh^- 



B.l. Locations 

The locations of the acoustic peaks in the CMB power spectrum depend on the value of 
sound horizon at decoupling epoch rs{7]dec) = Jq'''"" drj' Cs and the angular diameter distance to 
the last scattering surface, dA{zdec)- Comparing with numerical calculations it was shown (see 
(Efstathiou & Bond 1999; Hu et al. 2001; Doran & Lilley 2001) and references therein) that the 
spherical harmonic which corresponds to the m-th acoustic peak is well approximated by the relation 



_ dA{Zde c) 
rs{zde 



^Prn = {m-^m)^'^^^, (B2) 



where 0„i take into account the shift of m-th peak from its location in the idealized model which 
is caused by driving effects from the decay of the gravitational potential. Doran and Lilley (2001) 
give an accurate analytic approximation in the form 

^m = 4>- S(j)m , (B3) 

where 4> is overall phase shift of the spectrum (or the first peak) and S(j)m is a relative shift of each 
peak and dip caused by the Doppler shift of the oscillating fluid. For the overall phase shift of the 
spectrum they find 

4> = (1.466 - 0A66ns)air^^, (B4) 



where 

^ ^ _ 0.0416 f l + pjp^ \ f n Y f Zdec \ 

r*-Prad{Zdec)/Pm{Zdec)- ^ ^ gg^g J ^2. 726 J llOOoJ 

is the ratio of radiation to matter at decoupling, and 

ai = 0.286 + 0.626u;b , a2 = 0.1786 - 6.308^,, + 174.9u;^ - 1168u;^ 
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are fitting coefiicients. Here and below the numbers in the expressions are obtained for a present 
CMB temperature of Tq = 2.726K and the ratio of densities of massless neutrinos and photons 
Pi^/p-y = 0.6813 for three massless neutrino species (correspondingly fi, = Pu/ip-y + Pu) = 0.405). 
All values can be easily scaled to other values of Tq and f,^. 

The relative shift of the 1st acoustic peak is zero, (5^i = 0. For the 2nd one it is 

6(P2 = CO - cir* - C2/rT + 0.05(n, - 1) , (B5) 

with 

Co = -0.1+0.213e"^^'^^ ci = 0.015+0.063e-^^°°'^6 , C2 = 6-10-^+0.137(a;6-0.07)^ C3 = O.S+TOwft, 
and for the 3rd peak 

5^3 = 10 - di4^ + Om{ns - 1) , (B6) 

with 

di = 9.97 + 3.3iOb, d2 = 0.0016 + 0.196a;6 + 2.25 • lO'^w^^ 

The formula (B2) is correct also for the location of dips if we set m = 3/2 for the 1st dip and 
m = 5/2 for the 2nd dip. The relative shift of the first dip given by (Dor an &; Lilley 2001) is 

S(t)3/2 = bo + 61 exp + 0.158(n, - 1) (B7) 

with 

5o = -0.086 - 2.22wfe - 140^^ , h = 0.39 - 18.1a;b + 440w^, 62 = -0.57 - 3.8exp(-2365w^) . 
The angular diameter distance to the last scattering surface is given by 

Q 

dA{zdec) = — — ==X(% - %ec) , (B8) 
-nOV l^'fcl 

where x(^) = ^) sin a; or sinhx for flat, closed or open models respectively, and 



7o ^/^rad{z + 1)4 + ^rn{z + if + ^K + ^k{z + if 

Since, the sound speed in the pre-recombination plasma is 



= c/V3(l + R) with R = 3pb/Ap^ = 30315(ro/2.726)^'u;fea (BIO) 
and scale factor is well approximated by 
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with 

4.16-10-5 /I +p^/p-^W To y _ Veq 



OJ. 



m 



1.6813 J\2.726j ' ' 2(^2-1)' 



the integral for sound horizon can be reduced to the analytic formula 

19.9 / To \\ Vl + Rdec + y/Rdec + Req 



The deviation of the acoustic extrcma locations calculated using formulae (B2-B12) from the 
values obtained by CMBfast code is < 3% for a sufficiently wide range of parameters. 



B.2. Amplitudes 

The amplitude of the 1st acoustic peak can be approximated by the following expression 



— 



pi 



1 



27r 



(B13) 



where 



A = 0.838^ = exp [ai + 020;^^^ + azUcdm + 04^6 + o-b^h + aQ^hi^cdm + ai'^k + a^.'^l + a^ins - 1)] 

(B14) 

and CiJ^ is given by (A9). We have re-determined the best-fit coefficients Oj using the values of the 
1st acoustic peak amplitudes from CMBfast for the grid of parameters given below. Their values 
are 



ai = 2.326, 02 = 8.906, 03 = -7.733, 04 = -115.6, = 35.66, 
as = -7.225, ay = 1.96, ag = -11.16, ag = 4.439. 



(B15) 



The deviations of this approximation from the numerical value obtained by CMBfast are < 5% 
within the range of cosmic parameters, 0.2 < < 1.2, < < 0-8, 0.015 < < 0.12, 
0.8 < ra, < 1.2 and 0.4 < /i < 1.0. 



To calculate the amplitudes of the 2nd and 3rd peaks we use the relations 



Ho 



and = 



-P3 \^P3 



given by Hu et al. (2001). This leads to the following amplitudes 
Ap^ = Apj^H2i^m,^b,ns), with H2 



0.925a;{^^^2.4"--i 



1 + (u;6/0.0164) 



1/5 ' 



Apg = Ap^H2{Qm,^b,ns), with H3 



2.17(j^593.6"»-i 



1 + (a;b/0.044)^ [1 + 1.63(1 - a;b/0.071)a;^ 



(B16) 



(B17) 
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This approximation for Ap^ deviates by less than 5% from the value obtained with CMBfast 
for parameters within the range specified above. The accuracy of Ap^ is better than 9%. For some 
parameter values the second peak is under estimated leading to this somewhat poorer accuracy. 
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